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ABSTRACT 

Electron temperatures derived from the He I recombination hne ratios, des- 
ignated Te(He I), are presented for 48 planetary nebulae (PNe). We study the 
effect that temperature fluctuations inside nebulae have on the Te(He I) value. 
We show that a comparison between Te(He I) and the electron temperature de- 
rived from the Balmer jump of the H I recombination spectrum, designated 
Te(H I), provides an opportunity to discriminate between the paradigms of a 
chemically homogeneous plasma with temperature and density variations, and a 
two-abundance nebular model with hydrogen-deficient material embedded in dif- 
fuse gas of a "normal" chemical composition (i.e. ~ solar), as the possible causes 
of the dichotomy between the abundances that are deduced from collisionally 
excited lines to those deduced from recombination lines. We find that Te(He I) 
values are significantly lower than Te(H I) values, with an average difference of 
(Tc(H I)-Te(He I)) = 4000 K. The result is consistent with the expectation of the 
two-abundance nebular model but is opposite to the prediction of the scenarios 
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of temperature fluctuations and/or density inhomogeneities. Prom the observed 
difference between re(He I) and Te(H I), we estimate that the fiUing factor of 
hydrogen-deficient components has a typical value of 10~^. In spite of its small 
mass, the existence of hydrogen-deficient inclusions may potentially have a pro- 
found effect in enhancing the intensities of He I recombination lines and thereby 
lead to apparently overestimated helium abundances for PNe. 

Subject headings: ISM: general — planetary nebulae 

1. Introduction 

Two long-standing problems in nebular astrophysics are a) heavy element abundances 
relative to hydrogen obtained from coUisionally excited lines (CELs) are generally lower than 
those determined from optical recombination lines (ORLs); and b) electron temperatures de- 
duced from the coUisionally excited [O III] nebular-to- auroral forbidden fine ratio - hereafter 
Te{[0 III]) - are systematically higher than those determined from the Balmer jump (BJ) of 
H I recombination spectrum - hereafter Te(lH) - (see Liu 2003 for a recent review). 

Several solutions have been proposed for the discrepancies. Peimbert (1967) showed that 
temperature fluctuations within nebulae can lead to a Te{[0 III]) which overestimates the 
average emission temperature of the H I recombination spectrum and the [O III] nebular lines. 
As a consequence, the 0^"'"/H"'" ionic abundance ratio derived from the intensity ratio of the 
[O III] nebular lines to H/9, assuming Tc([0 III]) as the electron temperature, may be grossly 
underestimated. The presence of temperature fluctuations will also lead to a Te([0 III]) that 
is systematically higher than Te(H I), the temperature derived from the H I recombination 
spectrum. Evidence in favor of this was flrst presented by Peimbert (1971) who measured 
values of Te(R I) in several H ll regions and planetary nebulae (PNe) from the observed 
magnitude of the H I Balmer jump and found that they are systematically lower than the 
corresponding values of T'e([0 III]) of the same nebulae. A parameter t"^ was introduced by 
Peimbert (1967) to characterize the magnitude of temperature fluctuations in a given nebula. 
The scenario has however been found to be unable to account for the low abundances derived 
from temperature-insensitive IR CELs. High quality measurements of IR fine-structure lines 
with the Infrared Space Observatory (ISO) for a number of PNe have shown that these lines 
of very low excitation energies generally yield abundances similar to the values inferred from 
optical/UV CELs. Again they are lower with respect to ORL abundances (e.g. Liu et al. 
2000, 2004b; Tsamis et al. 2004). Moreover, recent observational studies of several samples of 
PNe show that typical values of t"^ derived by comparing T'e([0 III]) and Te(li I) are so large 
that they are well beyond the predictions of any photoionization models (Liu & Danziger 
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1993; Zhang et al. 2004) unless extra energy input (such as from magnetic reconnection or 
from shocks) is considered. 

A more plausible solution for the discrepancies of temperature and abundance determi- 
nations using CELs and ORLs is the two-component nebular model with hydrogen-deficient 
inclusions embedded in the diffuse nebula first proposed by Liu et al. (2000). Unlike 
the purely phenomenological temperature fluctuations scenario, the two-component nebu- 
lar model provides a physical explanation for the temperature and abundance determination 
discrepancies. The model assumes that there is a small amount of H-deficient material em- 
bedded in the diffuse nebula of "normal" composition. Due to the much enhanced cooling 
rates by IR fine-structure lines of heavy element ions as a result of the high metallicity, 
hydrogen-deficient inclusions have so low an electron temperature that they emit copiously 
in ORLs but essentially not at all in CELs. In this scenario, the high heavy elemental 
abundances yielded by ORLs simply reflect the enhanced metallicity in those H-deficient 
inclusions, rather than representing the overall chemical composition of the whole nebula. 
Detailed two-abundance photoionization models incorporating H-deficient inclusions have 
been constructed by Pequignot et al. (2003). The models provide a satisfying explanation 
for the apparent large discrepancies between the results of temperature and abundance de- 
terminations using CELs and ORLs, for the PNe NGC6153, Ml-42 and M2-36. Contrary 
to the temperature fluctuations scenario, the two-component nebular model does not need 
any extra energy input to reproduce the large temperature inhomogeneities within PNe. On 
the other hand, the presence of such H-deficient inclusions is not predicted by any of the 
current theories of the nucleosynthesis and evolution of low- and intermediate-mass stars 
and their nature and origin remain unclear. One possibility is that they are evaporating icy 
planetesimals (Liu 2003). 

More plasma diagnostic methods are helpful to probe nebular physical structure and 
to advance our understanding of the aforementioned fundamental problems. Earlier studies 
(e.g. Peimbert et al. 2000, 2002) show that temperature fiuctuations can affect the intensities 
of He I lines. Peimbert et al. (1995) presented a method to determine the average emission 
temperature of He I lines, - hereafter Tc(He I) - by adjusting the optical depth, the 
parameter and electron temperature to reconcile the He"'"/H"'" ratios yielded by individual 
He 1 lines. The method however requires accurate atomic data in addition to high quality 
line intensity measurements. 

Increasingly reliable atomic data are becoming available. A calculation for the effective 
recombination coefficients of a selection of He 1 lines was given by Smits (1996). Sawey & 
Berrington (1993) studied the contributions of coUisional excitation from the He° 2s^S and 
2s^S meta-stable levels by electron impacts to the observed fluxes of He I lines. Combining 
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the recombination data of Smits (1996) and the collision strengths of Sawey & Berrington 
(1993), Benjamin et al. (1999) presented improved He I line emission coefficients and fit the 
results with analytical formulae. More recently, Benjamin et al. (2002) studied the effects 
of the optical depth of the 2s"^S meta-stable level on He I line intensities. The availability of 
these improved atomic data has made it possible to obtain secure measurements of Te(He I), 
given high quality spectroscopic data. 

The main goal of the current paper is to demonstrate that He I emission lines provide a 
valuable probe of nebular thermal structure and a method to discriminate different scenarios 
proposed as the causes of the CEL versus ORL temperature and abundance determination 
discrepancies. In Section 2, we present electron temperatures derived from He I lines and 
confront the results with the predictions of the scenarios of temperature fluctuations, density 
inhomogeneities and of the two-abundance nebular model in Section 3. We summarize our 
results in Section 4. 



2. Electron temperatures from He I lines 

2.1. Method and atomic data 

Benjamin et al. (1999) provide analytic formulae for the emissivities of He I lines as 
a function of electron temperature for = 10^, 10"^ and 10^ cm~^, including contributions 
from recombinations of He^ with electrons (Smits 1996) and from the effects of collisional 
excitation from the 2s^S and 2s^S meta-stable levels (Sawey & Berrington 1993). Using their 
formulae, the intensity ratio of two He I lines is given by 

£l = ^lr,'r'^exp(^lp), (1) 

where Tc4 = Tg/lO'^ K. Values of the fitting parameters Oj, bi and Cj are given by Benjamin et 
al. (1999) for the temperature range 5 000-20 000 K, within which the emissivity fits have a 
maximum error of less than 5%. For lower temperatures, we have obtained similar fits using 
the numerical values given by Smits (1996) and Sawey & Berrington (1993) and present the 
results in Appendix A. For singlet lines. Case B recombination was assumed. 

By measuring the intensity ratio of two He I recombination lines, one can in principle 
determine the electron temperature from Eq. (1). This is however not an easy task given 
the relatively weak dependence of the line ratio on temperature. Apart from recombination 
and collision excitation, other minor effects that can potentially affect the observed line 
intensities may also need to be considered, such as self- absorption from the meta-stable 2s ^S 
level which can be significantly populated under typical nebular conditions. In order to 
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achieve the best measurements of the average He I hne emission temperature, one needs to 
take into account a variety of considerations, including the strengths and abihty to measure 
the hues involved in the ratio, the accuracy and reliability of the atomic data for the lines, 
the sensitivity of the ratio to Tg, and the magnitudes of contaminations by other unwanted 
side effects such as collisional excitation and self-absorption. After a thorough investigation, 
we decided that the best line ratio suitable for temperature determinations is probably the 
A7281/A6678, for the following reasons: 1) The AA7281, 6678 hues are amongst the strongest 
He I hues observable in the optical, after the A5876 and A4471 hues, and both hues fall 
in a spectral region clear of serious telluric absorption/emission and blending by any other 
known strong spectral features (see Fig. 1 for a typical nebular spectrum showing the relevant 
spectral region); 2) The atomic data associated with the AA7281, 6678 lines arc amongst the 
best determined; 3) The He I A7281/A6678 ratio is more sensitive to temperature but less 
sensitive to density compared to other ratios of strong He I lines; 4) Both the A6678 and 
A7281 lines arc from singlet states and thus arc essentially unaffected by the optical depth 
effects of the 2s '^S level; 5) Given the small wavelength span between the AA7281, 6678 lines, 
measurements of their intensity ratio are less sensitive to any uncertainties in reddening 
corrections and flux cahbration (c.f. also Liu et al. 2004b). 



2.2. Results 

In this work, a total of 48 PNe are analyzed. Amongst them, 23 were observed with the 
4.2 m William Herschel Telescope at La Palma and the ESO 1.52 m telescope, as described 
in Zhang et al. (2004). For the remaining sources, spectral data are from Liu et al. (2004a), 
Tsamis et al. (2003a), Wesson et al. (2004), Ruiz et al. (2003) and Peimbert et al. (2004). 
The corresponding references are listed in the last column of Table 1. 

In Fig. 2, we plot the He I A7281/A6678 ratio versus re(He I A7281/A6678) for different 
electron densities. The observed He I A7281/A6678 ratios along with measurement uncer- 
tainties are also overplotted. The electron densities of individual nebulae are taken from 
Zhang et al. (2004), Liu et al. (2004a), Tsamis et al. (2003a), Wesson et al. (2004), Ruiz 
ct al. (2003) and Peimbert et al. (2004). The observed line fluxes have been corrected for 
dust extinction using reddening constants taken from the same references. The resultant 
electron temperatures derived from the A7281/A6678 ratio are presented in Table 1. The 
associated uncertainties were estimated based on the line ratio measurement errors. Inspec- 
tion of Fig. 2 shows that most PNe in our sample have Te(He I) less than 10 000 K. For this 
low temperature regime, the A7281/A6678 ratio is quite insensitive to electron density. Thus 
our temperature determinations should be hardly affected by any uncertainties in density 
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measurements. 

For comparison, we have also determined temperatures using the He I A7281/A5876 
ratio. For this ratio, additional uncertainties may arise due to the optical depth effects of 
the 2s level on the intensity of the A5876 line. Based on the hne emissivities given by 
Benjamin et al. (2002), we have also plotted in Fig. 3 one curve showing the variations of 
the A7281/A5876 ratio as a function of for the case of A^'c = lO^cm"^ and rsggg = 100, 
the optical depth at the centre of the He I 2s^S - 3p^P° A3889 line. Fig. 3 shows that 
uncertainties in r^^sg may in principle lead to an error of approximately 1000 K in temperature 
determinations and that the uncertainty increases with increasing temperature. Values of 
Te(He I) derived from the A7281/A5876 ratio are also tabulated in Table. 1. Possible errors 
caused by uncertainties in T3889 have not been included in the error estimates. 

In Fig. 4, we compare electron temperatures derived from the A7281/A6678 ratio and 
from the A7281/A5876 ratio. Inspection of the figure shows that temperatures derived from 
the two ratios are generally in good agreement. There is however some evidence that the 
A7281/A5876 ratio tends to yield lower temperatures. Two effects may be responsible for the 
small offsets. In our calculations, the possible effects of self-absorption on the intensity of the 
A5876 line were not considered. If T3889 > 0, then our calculations may have systematically 
underestimated the A7281/A5876 temperatures (c.f. Fig. 3). Alternatively, if there is some 
small departure from pure Case B to Case A recombination for the He I singlet lines, then 
our calculations may have underestimated temperatures derived from both ratios. But this 
efi^cct is expected to be small for temperatures derived from the A7281/A6678 ratio as both 
lines involved in the ratio are from singlet spectral terms, and as a result the effect may have 
largely canceled out. It is difficult to discriminate between these two possibilities. In any 
case, the relatively small offset between the two temperatures as shown in Fig. 4 suggests 
that both effects are not large and should not have affected our temperature determinations 
by a significant amount. 

Table 1 also fists values of Te(H I), the temperature derived from the ratio of the 
hydrogen Balmer discontinuity at 3646 A to Hll A3770, taken from Zhang et al. (2004). A 
comparison between Te(He I) and Te(H I) is given in Fig. 5. It shows that, with the only 
exception of NGC40, Te(He I) is lower than Tc(H I). The average difference between the two 
temperatures for the whole sample is (re(H I)-Te(He I))=4000K. 

A measure of the electron temperature characterizing the He I zones can also be obtained 
by observing the He I discontinuity at 3421 A, produced by He+ recombinations to the He I 
2p3po ig^gi (^Lj^ ^ Danziger 1993). With re(H I) and Tc(He I) as free parameters, the H I 
and He I recombination continuum spectrum is calculated to fit the observed continuum 
spectrum of NGC6153. The best fit yields Te(He I) = 4000 ± 1500 K (Fig. 6), consistent 
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within the uncertainties with the value 3350± 1000 K derived from the He I A7281/ A6678 hne 
ratio. The He I A3421 discontinuity is much weaker than the H I Balmer discontinuity and 
falls in a spectral region of even shorter wavelengths crowded by strong emission lines such 
as the O III Bowen fluorescence lines AA3428,3444, and is therefore much more difficult to 
measure. Nevertheless, high quality spectroscopic observations of this discontinuity should 
be invaluable. 



3. Discussion 

3.1. Temperature fluctuations and density inhomogeneities 

The discrepancy between electron temperatures derived from the [O III] forbidden line 
ratio and those deduced from the hydrogen Balmer jump was first discovered by Peimbert 
(1971) and was ascribed to temperature fiuctuations within nebulae. For a nebula with 
small magnitude temperature fiuctuations, for a given ionic species of number density Ni, 
the nebular thermal structure can be characterized by an average temperature Tq and a 
mean square temperature fiuctuation parameter defined as (Peimbert 1967), 



and 



.2,.M _ !{T,-TofN,N4V 



For Tc(H I) derived from the H I recombination spectrum Balmer discontinuity, it can 
be shown that (Peimbert 1967), 

Te(HI) =To(l-1.67i'). (4) 

The deduction implicitly assumes that <C 1. 

Adopting the analytic expression for the fit of line emissivity as described in Section 2.1, 
we can write the fiux of a He I recombination line as, 

/(He I, A,) = / N,N{)le+)aiTll e^^{^)dV. (5) 



Then if ^ 1, one can rewrite the expression as a Taylor series expanded around the 
average temperature Tq, keeping only terms to the second order (Peimbert 1967), 
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where T04 = Tq/IO'^K. Then from Eqs. (2), (3) and (6), we have, 

/(HeI,Ai) _ oi T \ /\ \b1-b2 r C1-C2 

T \ ^ = — Te4 He I, A1/A2K exp 
/(He I, A2) a2 



T^{Re I, A1/A2) 

0'2 V -^04 / 



where 



^(^o) = \ {61(61 - 1) - 62(62 - 1) - 2[ci(6i - 1) - C2(62 - 1)]^ + ^i^} 
Introducing a parameter B and when <^ 1, we can obtain from Eqs. (7) 



Te4(He I, Ai/A2)^-^^ exp ^ [T,,{1 - Bef'-'^ exp 



Te4(He I,. 

B = ro4^(ro)/[ci - C2 - (61 - 62)704]. 

Thus if < 1 and l^t^l < 1, we have 
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(7) 



(8) 



(9) 



(10) 



In Eq. (10) the coefficient in front of is a slowly varying function of Tq. For Tc(He I; 
A1/A2) derived from the A7281/A6678 ratio, according to the He I line emissivity fit param- 
eters 61, 62, Ci and C2 given by Benjamin et al. (1999) for between 5000 and 20 000 K 
and A^e = 10^ cm~^, this coefficient varies between 0.97 and 1.25. For < 5000 K and 
A^e = 10'^ cm~^, using the fit parameters given in Appendix A (Table 2), we find that this 
coefficient varies between 1.03 and 1.22. Therefore, as an approximation, we adopt a value 
of 1.07 for this coefficient, as calculated for = 10000 K and A^e — 10^ cm~^, and have 



re(He I, A7281/A6678) « To(l - imt^). 



(11) 



Under the scenario of small temperature fluctuations {t^ <^ 1), we can reasonably assume 
that To(He I) 7o(H I) = Tq (c.f. Section 3.2 for further discussion on the legitimacy of this 
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assumption) and i^(He I) fa i^(H I) — t"^. Consequently, a comparison of Eqs. (4) and (11) 
shows that Tc(He I, A7281/A6678) > Te(H I), which is exactly opposite to what is observed 
(c.f. Section 2.2, Table 1). 

Using Eqs. (4) and (10), we plot in Fig. 5 re(H I) as a function of re(He I, A7281/A6678), 
for the case of t'^ — 0.00, 0.04, 0.10 and 0.16, assuming a constant density of lO^cm"^. 
Obviously, the differences between Tc(H I) and Te(He I) are in the opposite direction to 
what is predicted by temperature fluctuations. It follows that the He I lines may arise 
from regions characterized by significantly lower temperature than that of the H I emission 
regions, i.e., To(He l)< To(H I). This, however, will invalidate the condition ^ 1, and thus 
is beyond the scope of the description of the temperature fiuctuations as originally envisioned 
by Peimbert (1967, 1971). A detailed discussion will be given below in Section 3.2. 

Another theory that has been proposed to explain the ORL/CEL abundance determi- 
nation and the BJ/CEL temperature determination discrepancies is density variations in a 
chemically homogeneous nebula (Rubin 1989; Viegas & Clegg 1994). In this scenario, be- 
cause the [O III] A4363 auroral line has a much higher critical density than the AA4959,5007 
nebular lines and is therefore less affected by coUisional de-excitation in high density re- 
gions, the presence of such high density regions will lead to an overestimated T'e([0 III]) and 
consequently underestimated CEL abundances, analogous to the effects of temperature fluc- 
tuations. The presence of moderate density inhomogeneities in PNe is confirmed by recent 
observations (e.g. Zhang et al. 2004). However, given that H I and He I recombination lines 
are permitted transitions and therefore hardly affected by collisional de-excitation, such a 
scenario predicts Te(He I)~ Te(H I) (as long as He I zone ~ H I zone; see discussion below in 
section 3.2), also inconsistent with observations. We thus conclude that for a chemically ho- 
mogeneous nebula, the possible presence of large density inhomogeneities also fails to explain 
the observed differences between Te(lil) and T'e(He I). 



3.2. Two-abundance nebuletr models 

The two-abundance nebular model proposed by Liu et al. (2000) assumes the existence 
of another component of H-deficient, ultra-cold (Te ~ 10^ K), ionized gas embedded in 
the diffuse gaseous nebula of "normal" (i.e. ~ solar) chemical composition. In this model 
emission from the H-deficient ultra-cold ionized regions has a larger contribution to the He I 
than to the H I lines, which are still dominated by emission from the normal component 
under a typical temperature of Tc ~ 10^ K. Thus, while the model predicts that Te(H I) 
< Tc([0 III]), it also predicts that Te(He I) < Te(H I), in agreement with observations (Liu 
2003; Fig. 5). 
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Quantitatively, in a two-abundance nebula model, the intensity ratio of two He I lines 
can be written as, 



7(He I, Ai) ^ [N,N{Re-^)eiV]h + [N,N{}ie-^)e,V]i 
/(He I, A2) [N,N{Re+)e2V]h + [N,N{Re+)e2V]i ' 



(12) 



where £1 and 82 are the He I line emission coefficients, V is the volume of the emitting 
regions, and I and h refer to the low-metallicity regions (i.e. the diffuse nebula of "normal" 
composition) and the high-metaUicity regions (i.e. the cold H-deficient inclusions). Then 
using the analytic fit to the He I line emissivity, we have. 



7(HeI,Ai) _ ai 



/(He I,A2) 



= ^re4(HeI,Ai/A2)''i-''^exp 



02 



Cl - C2 



fli^ei exp ( 



Te4(He I,Ai/A2 



ai^ci cxp ( 



a2T^l exp (^)] ^ + [a2T^! exp 



(13) 



where /ig = {Ng)^/ {Ne)i and /^He = [-^(He^)]/i/[A^(He "*")]; are the electron density and He"*" 
ionic density contrasts between the two nebular components, respectively, and uj — Vh/Vi 
is the volume filling factor of the H-deficient component. Considering that the component 
of ultra-cold inclusions are hydrogen-deficient, the H I recombination line and continuum 
emission is characterized by the high temperature of the normal component, i.e. Te(li l)f» 

Given (Te)^, //e, A*He a-nd cu, one can thus obtain a relation between Te(H I) and Te(He I). 
The current available observations are however not sufficient to provide full constraints of all 
four free parameters. Detailed photoionization models incorporating hydrogen-deficient in- 
clusions have been constructed by Pequignot et al. (2003) to account for the multi- waveband 
observations of PNe NGC6153, M 1-42 and M2-36. If we make the simplified yet not un- 
reasonable assumption that the postulated hydrogen-deficient inclusions in all PNe are of 
similar characteristics, then guided by the photoionization models of Pequignot et al. (2003), 
we assume that for all PNe analyzed in the current work, {Te)h — 1000 K, pie — 100 and 
//He — 25. Under these assumptions, Te(R I) as a function of Te(He I) is plotted in Fig. 5 for 
filling factors of a; = 0, 10~^, 10~^ and 10~^. Here line emissivity parameters (oj, 6^, Cj) for a 
density of A'e = 10^ and 10^ cm~^ have been adopted for the "Z" and "/i" components, respec- 
tively. Adopting different densities however hardly affects our results since the emissivities 
of He I recombination lines are almost independent of electron density. 

Fig. 5 shows that a small amount of hydrogen-deficient material with a filling factor 
CO ~ 10~^ can account for the observed difference between re(He I) and Te(H I) for most 



- 11 - 



of the sample PNe. Such a small amount of hydrogen-deficient material would be difficult 
to detect by direct imaging, especially in the strong background emission of the "normal 
component". Note that in the two-abundance nebular model, heavy element ORLs arise 
almost entirely from the cold H-deficient inclusions, and their average emission temper- 
ature, Te(CNONe ORLs) should represent the true electron temperature prevailing in this 
H-deficient component. Thus one expects Te(CNONe ORLs) < Te(He I) (Liu 2003), a predic- 
tion that has also been confirmed by the available limited measurements of Te(CNONe ORLs) 
(Liu 2003; Wesson et al. 2003; Tsamis et al. 2004; Liu et al. 2004b; Wesson et al. 2004). Fur- 
ther accurate measurements of re(CNONeORLs) are thus essential to better constrain the 
physical conditions prevailing in these regions. For this purpose, high S/N ratio, high spec- 
tral resolution (and preferably also high spatial resolution) spectroscopy on a large telescope 
is required. 

An important consideration is whether the presence of a substantial He^"^ zone in some 
PNe may contribute significantly to the observed large difference between Te(H I) and 
Te(He I). Due to a higher heating rate and less efficient cooling, the electron temperature in 
the He^^ zone is generally higher than that in the He^ zone. This is indeed the case found 
from photoionization models for two PNe in our sample (NGC 7662 and NGC 3918) (Harring- 
ton et al. 1982; Clegg et al. 1987). It follows, then, that the average temperature of the H+ 
zone is higher than that of the He"*" zone since the former encompasses both the He"*" and He^"*" 
zones. However, none of the photoionization models published so far that assume a chemi- 
cally homogeneous medium predict a huge Te(li I) — Te(He I) difference as much as 4000 K 
found in the current work. As an example, the classic photoionization model of NGC 7662 
constructed by Harrington et al. (1982) yields only minute values of t^(H I) and t^(He I) 
and, as expected, in this model, To(H t)^ Tc(H I) and To(He l)^ Tc(He I). The model yields 
To(H I) of 12 590 K, in good agreement with the value Tc(H I) of 12 200 ±600 K deduced from 
the H I Balmcr discontinuity. However, their model predicts a value of 11 620 K for To(He I), 
significantly higher than the value 7690 ± 1650 K derived from the He I A7281/A6678 line 
ratio. Finally, we note that the observed phenomenon of large difference between Te(H I) 
and T'e(He I) persists in a number of low excitation PNe in our sample, such as H 1-35, Hu 2-1 
and M 1-20 where the He^+ zone is essentially absent. As a consequence, we conclude that 
the presence of a Hc^+ zone in PNe of high excitation is unlikely to play a major role in 
causing the large Te(H I), re(He I) discrepancy. 

In Fig. 7, we plot T'e(H I) and Te(He I) against the excitation class, E.G., calculated 
following the formahsm of Dopita & Meatheringham (1990) and tabulated in Table 1. Fig. 7 
shows that there is a positive correlation between Te(H I) and the E.G.. Excluding two 
peculiar PNe, He 2-118 and M2-24, we obtain a linear correlation coefficient of 0.66. This 
is consistent with the expectation that a hotter central star produces a higher heating rate 
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per photoionization, resulting in a higher nebular temperature. Fig.7 also shows a weaker 
correlation between Te(He I) and the E.G.. Again after excluding NGC40 in addition to 
He 2-118 and M2-24, the data yield a correlation coefficient of 0.41 between Te(He I) and 
the E.G.. Note that in the two- abundance nebular model, He I lines are strongly enhanced 
by emission from the H-deficient inclusions. As a consequence, the E.G. is no longer a 
good indicator of Tc(He I). Rather, in this scenario, Te(He I) will be mainly determined by 
the amount and properties of the postulated H-deficient inclusions, including their spatial 
distribution, density, chemical composition and filling factor. 

It is interesting to note that amongst the sample PNe the only nebula that exhibits 
a higher value of Te(He I) compared to Tc(H I), i.e. opposite the predictions of the two- 
abundance model, is NGG40 where re(He I) is about 3500 K higherth&n T^(R I) (c.f. Fig.5). 
NGG40 is ionized by a WG8 Wolf-Rayet central star (Smith 1969; Growther et al. 1970). It 
is possible that in this particular object the nebular emission line spectrum is contaminated 
by emission from the strong stellar winds from the high-temperature H-deficient envelope of 
the central star, leading to apparently higher Te(He I) with respect to Tc(H I) (c.f. Liu et al. 
2004b). Further observations of this peculiar PN are needed to clarify the situation. 



All He I hnes observable in the optical and UV are essentially entirely excited by radia- 
tive recombination. As a consequence, their intensities can be significantly enhanced by the 
presence of a small amount of ultra-cold, H-deficient inclusions, leading to an overestimated 
overall He/H abundance characteristic of the entire ionized region (Liu 2003, 2004). For a 
large sample of PNe, Zhang et al. (2004) showed that there is a positive correlation between 
He/H abundance and the difference between the temperature derived from the [O III] for- 
bidden lines and from the hydrogen Balmer discontinuity, lending further support to this 
possibility. In the following, we present an analytical method to obtain a quantitative esti- 
mate of the possible enhancement of the He/H abundance due to the postulated presence 
of H-deficient inclusions in PNe. The assumptions about the two-abundance nebular model 
are the same as those described above in Section 3.2. 

The intensity ratio of a He I recombination hne to that of H/3, /(He I, Aj) /I (HP) observed 
from a nebula is given by 



3.3. 



Helium abundance 



/(He I, A,) 



Ar(He+,A,)(£,)z4861 

iV(H+)(£H/3)Ai 

[N^N{Re+)eiV]h+[N^N{Re+)eiV]i 4861 
[N,N{m)eupV]h + [N,N{m)enpV]i A, ' 



(14) 
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where A'"(He''", Ai)/A'"(H+) is the He+/H+ ionic abundance ratio determined assuming a chem- 
ically homogeneous nebula. Prom the above equation we find, 



Ar(He+) 
Ar(H+) 



(15) 



For the small amount of H-deficient material as hypothesized in the two-abundance 
model, the average He+/H+ ionic abundance ratio for the entire nebula is essentially identical 
to the value for the diffuse material, i.e. [A^(He^)/A^(H+)]/. In such a case, the traditional 
method of abundance analysis assuming a chemically homogeneous nebula will overestimate 
the He/H abundance by a factor of [co'/Xe/XHe(^i)/i/(^i)i + !]• For the representative values of 
He — 100, Hiie — 25 and u; — 10"^, the He+/H+ ionic abundance ratio will be overestimated 
by a factor of 1.25. For example, the empirical analysis of NGC6153 by Liu et al. (2000) 
assuming a chemically homogeneous nebula yields a He/H elemental abundance ratio of 
0.136, which is a factor of 1.35 higher than the average He/H abundance ratio for the 
entire ionized region derived from the detailed two-abundance photoionization modeling of 
this nebula by Pcquignot et al. (2003). The much lower He/H abundance ratio obtained 
by two-abundance photoionization modeling is supported by observations of the coUisional 
excitation dominated He I 2s^S - 2p^P° A10830 fine (Liu 2003). 

Many PNe are known to show enhancement of helium with respect to the Sun. For 
example, the average He/H ratio obtained by Kingsburgh & Barlow (1994) for a large sample 
of Galactic PNe is a factor of 1.35 solar. The enhancement is often ascribed to the second 
and third dredge-up processes that occur during the post-main-sequence evolution stages of 
the progenitor stars of PNe. However, the extremely high He abundances deduced for PNe 
such as He 2-111 (Kingsburgh & Barlow 1994) are difficult to explain by the current theories 
of nucleosynthesis and dredge-up processes. It is possible therefore that the very high He/H 
abundances observed in those extreme He-rich PNe are actually caused by the contribution 
of H-deficient inclusions embedded in the nebulae. 

The CEL/ORL abundance discrepancy and the BJ/CEL temperature discrepancy which 
are ubiquitously observed amongst PNe are also found in H II regions (e.g. Tsamis et al. 
2003b). The determination of Te(H I) for H II regions is much more difficult than for PNe due 
to the generally much lower surface brightness of H II regions and the strong contamination 
of scattered star light to the nebular spectrum in dusty H II regions. Recently, Garcia-Rojas 
et al. (2004) presented deep echelle spectroscopy of the Galactic H II region NGC3576, 
covering the wavelength range from 3100-10400 A. They derived re(H I) of 6650 ± 750 K 
and T'e([0 III]) of 8500 ± 50 K, in agreement with the general pattern observed amongst 
PNe that Te(li I) is systematically lower than Te{[0 III]). On the other hand, from the He I 
AA6678, 7281 fine fiuxes reported in their paper, we obtain re(He I) of 6800 ± 600 K, which 
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is consistent with Te(H I) within the errors. Note that NGC3576 has a relatively small 
ORL/CEL abundance discrepancy (a factor of 1.8; Tsamis et al. 2003a; Garcia-Rojas et al. 
2004). Measurements of Te(H I) and Te(He I) in more H ll regions are needed, especially for 
those exhibiting large ORL/CEL abundance discrepancies in order to have a better picture 
of this problem in H II regions. 

For a few metal poor cxtragalactic H II regions, Peimbert et al. (2002) find that the 
electron temperatures derived from He I recombination lines are also systematically lower 
than the corresponding values deduced from the [O III] forbidden line ratio. If we ascribe 
the discrepancy to the presence of H-dcficicnt inclusions in these nebulae, then the effects of 
these inclusions on the determination of the primordial helium abundance can be estimated 
roughly from Eqs. (13) and (15). Using the line fluxes presented by Peimbert et al. (2002) 
and assuming that the H-deficient material in these metal-poor H ll regions have similar 
composition as those postulated to exist in PN NGC 6153, we estimate that the primordial 
He abundance could have been overestimated by as much as 5%. 

4. Conclusions 

In this paper, we show that He I recombination lines can be used to probe nebular 
thermal structures and provide vital information regarding the nature and physical causes of 
the long-standing ORL/GEL temperature and abundance determination discrepancies. We 
present electron temperatures derived from the He I A7281/A6678 ratio for 48 PNe. The 
He I temperatures are found to be systematically lower than those deduced from the Balmer 
jump of the H I recombination spectrum. The result is exactly opposite to the predictions of 
the scenarios of temperature fiuctuation and density inhomogeneities but in good agreement 
with the expectations of the two-abundance nebular model proposed by Liu et al. (2000). 
We estimate that a filling factor of the order of 10~^ of the ultra-cold, H-deficient material 
is sufficient to explain the observed differences between the H I and He I temperatures. We 
show that the possible presence of H-deficient inclusions in PNe may indicate that current 
estimates of He/H abundances in PNe may have been overestimated by a factor of ~ 1.25. 
The possible existence of hydrogen-deficient inclusions in H II regions could also cause the 
primordial helium abundance determined from metal-poor H II galaxies to be overestimated. 
We stress the importance of high S/N ratio, high resolution spectroscopy in the tackling of 
these fundamental problems. 

We thank the referee. Dr. V. Luridiana, for helpful suggestions and incisive comments 
that improved the paper significantly. YZ acknowledges the award of a Liu Yongling Scholar- 
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ship. Support for RHR was from the NASA Long-Term Space Astrophysics (LTSA) program. 

A. He I line emissivities (Te < 5000 K) 

Combining the He I recombination model of Smits (1996) and the coUisional excitation 
rates of Sawey & Berrington (1993), we have calculated emissivities for the He I AA4471, 5876, 
6678 and 7281 lines for < 5000 K - the temperature range not considered by Benjamin 
et al. (1999). At such low temperatures, the effects of coUisional excitation are however 
essentially negligible. We may write, following the analytic expression used by Benjamin et 
al. (1999), 

Si = ailll exp(cj/re4) ergcm^ s"\ (Al) 

where Te4 = Te/IO^K, and a^, hi and q are constants. The constants were derived using 
a least-squares algorithm. The results are tabulated in Tabic 2 for electron densities of 
A'^c = 10^, 10'^ and 10® cm^'^. The maximum errors of the fits for the temperature range 
312.5 — 5000 K, Cj, are also listed in the last column of the Table. 
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J I I I I I I I I I I I I I I L 

6000 6500 7000 

A (A) 



Fig. 1.— Spectrum of NGC7009 showing the He I AA5876, 6678 and 7281 hnes used to 
determine average emission temperatures of the He I recombination spectrum. The intensity 
was scaled such that -F(H/3) = 100. Note that Ha was heavily saturated in this spectrum. 
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2000 3000 4000 5000 7000 10000 13000 20000 



T^iHe I X7281/X6678)(K) 

Fig. 2. — The He I A7281/A6678 intensity ratio as a function of Te for electron densities from 
iVg = 10'^ to lO^cm"^. The observed ratios along with their uncertainties for our sample 
PNe are overplotted. 




Fig. 3. — The He I A7281/A5876 intensity ratio is plotted against for densities from 
Nc = 10^ to lO^cm"^. All curves are for zero A3889 line optical depth except one for which 
^3889 = 100 and = lO^cm"^ (see text for more details). 
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5000 10000 
T^{He I X7281/A6678)(K) 



Fig. 4. — Comparison of electron temperatures derived from tiie He I A7281/A6678 ratio 
with those derived from the He I A7281/A5876 ratio. The diagonal hne is a, y = x plot. 
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5000 10000 15000 20000 

7'e(He I X7281/X6678) (K) 

Fig. 5. — Tc(Hc I) versus Tc(H I). The solid diagonal line is a y = x plot. The dotted 
lines show the variations of Tc(H I) as a function of Tc(He I) for mean square temperature 
fluctuation parameter = 0.04, 0.10 and 0.16, assuming a constant electron density of 
10^ cm~^. The short-dashed, dot-dashed and long-dashed lines show the variations of Te(H I) 
as a function of re(He I) for u; — 10~^, 10~^ and 10~^, respectively, where u; is the volume 
filling factor of H-deficient inclusions. 
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Tg(H I)=6000 K 



I I I I I I I I I I I I 

3400 3600 3800 

A (A) 

Fig. 6.— Spectrum of NGC 6153 from 3250-3900 A, showing the He I discontinuity at 3421 A 
and the H I discontinuity at 3646 A. The dotted hne is a theoretical nebular continuum 
spectrum calculated assuming a He I temperature of 4000 K and a H I temperature of 
6000 K. For comparison, a theoretical spectrum assuming a He I temperature of 6000 K is 
also given (dashed hne), which clearly yields a poorer fit to the observed magnitude of the 
He I A3421 discontinuity than the dotted line. 
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Fig. 7. — a) Te(H I) and b) Te(He I) versus Excitation Class (E.G.). The solid lines are linear 
least-squares fits. 
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Table 1. Electron temperatures derived from He I line ratios. 



Source 


T,(He I A7281/A6678) 


T,,(He J A7281/A5876) 


T,,(H I) 


E.G. 


Ref. 


v^n z-i 


oyuu It iouu 


K/inn 4- 1 'inn 

04UU It ioUU 


1 nsnn -u finn 

lUoUU It DUU 


4, 


.87 




n i-oo 


RQnn 4-11 nn 
oyuu zt iiuu 


cccnn 4- 1 9nn 

OOUU It iZUU 


1 nnnn -i- i nnn 

lUUUU It IUUU 


2, 


.14 




ri i-ou 


conn 4- 1 9nn 
DoUU ± iZUU 


cc9nn 4- 1 Qnn 

OZUU it ioUU 


1 9i^nn 4- 1 nnn 

iZOUU it IUUU 


4, 


.99 




jne z-i±o 


7Qc;n 4-11 nn 


«ni^n 4-11 nn 

DUOU It iiUU 


1 SQnn 4- 1 nnn 

looUU It IUUU 


3 


.12 




nu i-z 


1 1 Knn 4- 1 ^nn 


1 1 7c:n 4- '^nnn 

11/ OU It OUUU 


9nnnn -i- 1 9nn 

ZUUUU It IZUU 


9, 


.69 


T C\A 


rlU Z- i 


OA nn 4-11 nn 

o4UU It i iUU 


7nnn 4- 1 9nn 

(UUU It iZUU 


1 nnnn 4- c;nn 
iuuuu It ouu 


1, 


.81 


VV U4 


TP 1 9Q7 


cnnn 4- i nnn 

OUUU It iUUU 


Konn 4- 1 nnn 

OZUU It iUUU 


1 nnnn -i- /inn 

lUUUU It 4UU 


6, 


.87 




TP 9nnQ 


Kfinn 4-11 nn 

OOUU It i iUU 


«7nn 4- 1 9nn 

u ( UU It iZUU 


1 1 nnn -u snn 

liUUU It oUU 


4, 


.78 


VV U4 


oODo 


oi nn 4- 1 nnn 

oiUU It iUUU 


7Qnn 4- T AKr\ 

1 oUU It i40U 


Qonn 4- nnn 
youu It yuu 


4, 


.59 


T n/1 

ijU4 


TP /1 1 Q1 


KKnn 4- 1 nnn 

OOUU It iUUU 


Knnn 4- i nnn 

OUUU It IUUU 


Qonn 4- 1 nnn 
yzuu It IUUU 


5, 


.12 


Tnq 
1 Uo 


TP A And 


Konn 4- 1 nnn 

OZUU It iUUU 


A 4_ 1 nnn 

400U It IUUU 


nqcn 4- 1 nnn 
yoou It IUUU 


4, 


.68 


Tnq 
1 Uo 


TP /IRQ/I 


r,A nn 4- 1 9nn 

04UU It iZUU 


A 1 nn 4- 1 nnn 

41UU It IUUU 


Qcnn 4- /inn 

oOUU It 4UU 


3 


.03 




TP /I77R 


«i en 4- 1 1 i^n 

DiOU It iiOU 


c;9nn 4- 1 Qnn 

OZUU It loUU 


Q«nn 4- A nn 

oDUU It 4UU 


3 


.08 




TP /IOQ7 


7Knn 4- 1 c;nn 
/ ouu It iouu 


«9Kn 4- 1 c;nn 

DZOU It IOUU 


1 n9nn -i- c;nn 
luzuu It ouu 


3 


.07 




TP K0^ 7 

OZi ( 


«i nn 4-11 nn 

DiUU It i iUU 


c;«nn 4- 1 9nn 

OuUU It IZUU 


1 9nnn 4- snn 

iZUUU It oUU 


5, 


.19 


VV U4 


i\/r 1 9n 

iVi i-ZU 


c;ssn 4- 1 9nn 

OooU It iZUU 


/I K c;n 4- 1 nnn 

400U It iUUU 


1 9nnn 4- 7nn 

IZUUU It (UU 


2, 


.97 




i\/r 1 /19 

iVi i-4Z 


oo«n 4- 1 nnn 

ZZOU It iUUU 


Kocn 4- 1 9nn 

OZoU It IZUU 


/innn -i- finn 

4UUU It OUU 


2, 


.21 




A/T 9 9/1 
iVi Z-Z4 


/( KKn 4- 1 c;nn 

400U It iOUU 


9Qnn 4- 1 nnn 

ZoUU It IUUU 


1 finnn 4- 9nnn 
lOUUU It zuuu 


1, 


.59 




iVi Z-OD 


97Qn 4- 1 nnn 
z / yu It IUUU 


Qi 9n 4- 1 nnn 

OlZU It IUUU 


«nnn -i- Ann 

OUUU It 4UU 


3 


.43 




i\/r Q 91 
iVi o-z± 


c;9nn 4- 1 9nn 

OZUU It iZUU 


QQ1 n 4- 1 nnn 
oyiu It IUUU 


1 1 nnn 4- /i nn 

IIUUU It 4UU 


A 

4, 


1 1 

.11 




i\/r Q Q9 

iVi O-OZ 


OA Qn 4- 1 nnn 

Z4oU It iUUU 


91 1 n 4- 1 nnn 

ZllU It IUUU 


A c;nn 4- i^nn 
40UU It OUU 


3 


.00 




Me 2-2 


6200 ± 1100 


6050 ± 1200 


11000 ± 500 


3 


.32 


W04 


NGC 40 


10240 ± 1900 


10580 ± 4200 


7000 ± 700 


0, 


.18 


L04 


NGC 3132 


8540 ± 1000 


8370 ± 1000 


10780 ± 1000 


3, 


.55 


T03 


NGC 3242 


4850 ± 1000 


4270 ± 1000 


10200 ± 1000 


5, 


.74 


T03 


NGC 3918 


6050 ± 1000 


5480 ± 1000 


12300 ± 1000 


6, 


.63 


T03 


NGC 5307 


6820 ± 1050 


6540 ± 1100 


10700 ± 1000 


5, 


.69 


R03 


NGC 5315 


5480 ± 1000 


5010 ± 1000 


8600 ± 1000 


3 


.57 


T03 


NGC 5315 


6310 ± 300 


5720 ± 400 


7500 ± 1000 


3 


.79 


P04 


NGC 5873 


5740 ± 1300 


5990 ± 1400 


12000 ± 500 


7, 


.02 




NGC 5882 


5340 ± 1000 


4660 ± 1000 


7800 ± 1000 


4, 


.73 


T03 
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Table 1 — Continued 



Source 


-Le[tie I A/ZOi/ADD/oj 


Je(,tie I A/zoi/Aoo/Dj 


^e(H I) 


E.G. 


Kei. 


NGC 


6153 


3350 ± 1000 


2700 ± 1000 


6000 ± 400 


3 


.83 




NGC 


6210 


6650 ± 1200 


6900 ± 1350 


8700 ± 800 


4, 


,69 


L04 


NGC 


6302 


7830 ± 1300 


6850 ± 1300 


16100 ± 500 


8, 


.06 




NGC 


6543 


6010 ± 1400 


5450 ± 1400 


6800 ± 400 


2, 


.94 




NGC 


6567 


7480 ± 1300 


6640 ± 1300 


12000 ± 800 


4, 


.12 




NGC 


6572 


8690 ± 1200 


7430 ± 1650 


10300 ± 1000 


5, 


.26 


L04 


NGC 


6620 


7660 ± 1300 


7120 ± 1300 


10000 ± 500 


5, 


.71 




NGC 


6720 


8120 ± 1700 


7810 ± 2200 


8000 ± 600 


5, 


.52 


L04 


NGC 


6741 


8500 ± 1300 


7340 ± 1880 


15000 ± 600 


6, 


.33 


L04 


NGC 


6790 


9910 ± 1500 


8790 ± 2000 


14000 ± 500 


6, 


.23 


L04 


NGC 


6803 


5000 ± 1100 


4500 ± 1200 


8500 ± 400 


4, 


.93 


W04 


NGC 


6833 


6000 ± 1100 


4000 ± 1200 


14000 ± 500 


3 


.38 


W04 


NGC 


6818 


3690 ± 1000 


3820 ± 1000 


12140 ± 1000 


7, 


.71 


T03 


NGC 


6818 


5310 ± 1300 


5880 ± 1300 


12500 ± 600 


7, 


.71 




NGC 


6826 


8290 ± 1500 


8520 ± 2000 


8700 ± 700 


3 


.26 


L04 


NGC 


6884 


7770 ± 1000 


6030 ± 1550 


11600 ± 400 


5, 


.34 


L04 


NGC 


7009 


5040 ± 800 


4620 ± 1100 


7200 ± 400 


5, 


.19 




NGC 


7027 


10360 ±1100 


9030 ± 2200 


12000 ± 400 


6, 


.97 




NGC 


7662 


7690 ± 1650 


7620 ± 1950 


12200 ± 600 


6, 


.79 


L04 



References. — [L04] Liu et al. (2004a); [P04] Peimbert et al. (2004); [R03] Ruiz et al. 
(2003); [T03] Tsamis et al. (2003a); [W04] Wesson et al. (2004) 
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Table 2. Fitting parameters for He I line emissivities. 



Emissivities 






bi 


Ci 


ei 


K-- 


= 10^ cm ^ 














^4471 


6.835x10" 


-26 


-0.8224 


-0.0074 


1.40% 




^5876 


1.838x10- 


-25 


-0.9745 


-0.0086 


1.11% 




J6678 


5.251x10- 


-26 


-0.9819 


-0.0088 


1.24% 




^7281 


9.104x10- 


-27 


-0.5594 


0.0007 


0.25% 


















J4471 


6.775x10- 


-26 


-0.8270 


-0.0041 


1.50% 




i5876 


1.824x10- 


-25 


-0.9702 


-0.0068 


1.14% 




J6678 


5.203x10- 


-26 


-0.9767 


-0.0068 


1.23% 




i7281 


9.463x10- 


-27 


-0.5289 


0.0068 


1.14% 




=10^ cm-3 














^4471 


6.746x10- 


-26 


-0.8337 


0.0067 


1.23% 




J5876 


1.819x10- 


-25 


-0.9509 


0.0035 


1.37% 




J6678 


5.185x10" 


-26 


-0.9581 


0.0035 


1.67% 




J7281 


9.663x10" 


-27 


-0.5307 


0.0179 


2.29% 
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